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Abstract 

A known generalization of the Stillinger-Lovett sum rule for a guest 
charge immersed in a two-dimensional one-component plasma (the sec- 
ond moment of the screening cloud around this guest charge) is more 
simply retrieved, just by using the BGY hierarchy for a mixture of sev- 
eral species; the zeroth moment of the excess density around a guest 
charge immersed in a two-component plasma is also obtained. The 
moments of the electric potential are related to the excess chemical 
potential of a guest charge; explicit results are obtained in several spe- 
cial cases. 
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1 Introduction 



One of us (L. S.) has derived a generalization of the Stillinger-Lovett sum rule 
for a guest charge immersed in a two-dimensional one-component plasma [1] : 
an exact simple expression for the second moment of the screening cloud 
around the guest charge was obtained, by using a mapping technique onto a 
discrete one-dimensional anticommuting-field theory. In the present paper, 
we first show that the same result can be obtained in a simpler way by just 
using the BGY hierarchy, which provides also more general results. 

The excess chemical potential of a guest charge (which can be expressed 
in terms of the charge density of the screening cloud) has an expansion in 
powers of the guest-particle charge Ze, which allows to compute the average 
of powers (moments) of the electric potential. 

We consider a classical (i.e. non-quantum) system of charged particles 
located in an infinite two-dimensional (2D) plane of points r e R 2 . According 
to the laws of 2D electrostatics, the particles can be thought of as infinitely 
long charged lines in the 3D which are perpendicular to the 2D plane. The 
electrostatic potential v at a point r, induced by a unit charge at the origin 
0, is thus given by the 2D Poisson equation 

Av(r) = -27r£(r). (1.1) 

The solution of this equation, subject to the boundary condition Vv(r) — ► 
as |r| — > oo, reads 

v(r) = -ln(£), (1.2) 

where r = |r| and the free length constant L, which determines the zero point 
of the potential, will be set for simplicity to unity. The Fourier component 
of this potential v(k) ccl/k 2 exhibits the characteristic singularity at A; = 0, 
which maintains many generic properties (like screening) of "real" 3D charged 
systems. 

A general Coulomb system consists of M mobile species a = 1, 2, . . . , M 
with the corresponding charges e a (which may be integer multiples of the 
elementary charge e). Mobile particles may be embedded in a fixed uniform 
background of charge density p b . The most studied models are the one- 
component plasma (OCP), which corresponds to M = 1 with e\ = e and 
Pb of opposite sign, and the symmetric two-component plasma (TCP), which 
corresponds to M = 2 with e\ = e, e 2 = — e, p b = 0. The interaction energy 
of a configuration {r-j, e ai } of the charged particles plus the background is 

E = Y^ e OLi e ajV (|r; - Tj\) +^e (Xi (j) b {r i ) + E b _ b , (1.3) 

i<j i 

where 0&(r) is the one-body potential created by the background and the 
background-background energy term E b -b does not depend on the particle 
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coordinates. In the case of point particles, for many- component systems with 
at least two oppositely species, the singularity of the Coulomb potential (jl.2p 
at the origin r = prevents, for small enough temperatures, the thermody- 
namic stability against the collapse of positive-negative pairs of charges. In 
those cases, one introduces to v a short-range repulsion which prevents the 
collapse. 

The Coulomb system is considered in thermodynamic equilibrium, at in- 
verse temperature (3 = l/{k&T). The thermal average over an infinite neutral 
system will be denoted by (• • •). In terms of the microscopic density of par- 
ticles of species a, n a (r) = J2i^a,ai^(r — r^), the microscopic total number 
density and the microscopic total charge density are defined, respectively, by 

^( r ) = H^*(r), p(r) = ^e a h a (r) + p b . (1.4) 

a a 

The microscopic electrostatic potential created by the particle-background 
system at point r is given by 

4>(r) = J dr'<u(r-r')p(r'). (1.5) 

At the one-particle level, the homogeneous number density of species a and 
the total particle number density are given respectively by 

n a = (n a (r)), n=(n(r)). (1.6) 

The charge density p = (p(r)) vanishes due to the charge neutrality of the 
system. At the two-particle level, one introduces the translationally invariant 
two-body densities 

n£,(|r-r'|) = 

= (n a (r)n a /(r')} - (n a (r))5 a>al 8(r - r'). (1.7) 
It is useful to consider also the pair distribution functions 

w(lr-r1) = n -' (|r ~ r ' D , (1.8) 

the (truncated) pair correlation functions h aa i = g aa > — 1, as well as the 
three-body analogous quantities 



, r yy/ ) = 4 3 L»(r,r',r") 

n a n a in a ii 



and the (truncated) three-body correlation function 

hal'ai'ir, r ' ; r ") = & a »{r, r', r") - h aa ,(\r - r'|) - h a , a// (\r' - r"\) 

-Va(|r"-r|)-l. (1.10) 
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The paper is organized as follows. In Section 2, we use the BGY hierarchy 
for studying the general mixture of M species of mobile particles embedded in 
a fixed uniform background. By taking the limit of one of the densities going 
to zero, we get the case of a guest charge. We retrieve the second moment 
of the screening cloud around a guest charge immersed in an OCP; we get 
also the zeroth moment of the excess total number density around a guest 
charge immersed in a TCP. In Section 3, the general formalism for relating 
the moments of the electric potential to the excess chemical potential of a 
guest charge is established. The following Sections study special cases when 
explicit calculations are possible: the high-temperature (Debye-Hiickel) limit 
in Section 4, the OCP at f3e 2 = 2 in Section 5, the TCP in Section 6. Section 
7 is a Conclusion. 

2 Sum rules for a guest charge immersed in 
a Coulomb system 

We wish to rederive and extend the result of [1] about the 2D OCP in which 
a point guest charge Ze is immersed. Let the charge density at r knowing 
that there is a guest charge Ze at the origin be p(r|Ze, 0). In [TJ, its second 
moment was shown to be 



J dr r 2 p(r\Ze, 0) 



/ r i — . / - 

7r/3 'en 



Z\l-^—\+Z 



(2.1) 



Our rederivation uses only the BGY hierarchy. 



2.1 General sum rule for a mixture with a background 

We start with the mixture of M mobile species, with a fixed uniform back- 
ground, described in the Introduction. Finally, we shall consider a mixture 
of only 2 species with respective charges e\ = e and e 2 = Ze\ at the end, the 
density n 2 will be chosen as 0, leaving only one guest charge. But for being 
able to consider the TCP as well, we start with the more general case of M 
mobile species. The neutrality constraint is 

J2n a e a = -p b . (2.2) 

a 

The mixture with a background has been studied in three dimensions by 
Suttorp and van Wonderen [2J. They used the BGY hierarchy and thermody- 
namical properties of the system for deriving, among other things, a second- 
moment sum rule, which however involves some thermodynamical functions 
(the partial derivatives of each density n a with respect to the background 
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density n&); there is no explicit expression for these partial derivatives. For- 
tunately, we found that, in two dimensions, the formalism becomes much 
simpler and only the BGY hierarchy has to be used (the thermodynamical 
properties are no longer involved). 

The second equation of the BGY hierarchy [3], with h aa >(r) the correlation 
function between a particle of species a at r and a particle of species a' at 
the origin, is 

/rV/w(r) = 

- J2 n °<" / dr " h a , a ii{r")e a e a »Vv{\r - r"|) 

- K a i(r)e a e a >Vv(r) - e a e a >Vv(r) 

J2n a >< fdr"h^ lafl (r,0,r")e a e al/ Vv(\r-r"\). (2.3) 

a" J 

The integral in the first term in the rhs of (12.31) is proportional to the electric 
field at r due to the charge distribution h a > a rr which has a circular symme- 
try around the origin. Thus, using Newton's theorem, one can rewrite this 
integral as 

/ dr" W(r")V<|r - r"|) = Vv(r) f dr" W(0- (2-4) 

J Jr"<r 

The integral in the rhs of (12. 4p can be written as J r » <r ■ ■ ■ = / . . . — f r '/ >r ■ ■ ■ 
and the perfect screening of the charge e a > gives [1] 

^e a iin a ii I dr" h a > a »(r") = -e a >. (2.5) 

a" 

Therefore (12.31) can be rewritten as 

/r x v/w(r) = 

e a n a tie a iiVv (r) / dr" h a i a " Or") 

Jr">r 

- K a i(r)e a e a <Vv(r) 

e a ^n Q » f dr"h® a ,,{r,0,r")e a ,,Vv{\r-r"\). (2.6) 

a" J 

In order to make a second moment to appear, we take the scalar product 
of both sides of (12. 6p with r and integrate on r. Integrating by parts the lhs 
and performing the integration on r first in the first term of the rhs, one finds 
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—2(3 1 J dr h aa i{r) = — ne a n a ne a n J dr r 2 h a i a n(r) 

a" 

+e a e a > j dr h aa r(r) (2.7) 
+e a n a "e a " J dr d(r" - r) ?4a'a"( r > °> r ") 



r - r") ■ r 



(r- r") 2 



(in the last term, since depends only on r and the difference v" — r, we 
have replaced the integration on v" by an integration on r" — r). An important 
simplification has occurred in 2D where r • Vt>(r) has the constant value — 1, 
while in three dimensions, with the potential v(r) = 1/r, one finds —v(r), a 
result which has led to a more complicated calculation in [2]. 

Now, we multiply both sides of (12. 7p by n a and sum on a. The term 
involving can be simplified by using symmetries under permutations of 
the variables. Indeed, has the symmetry property 

h^, a „(r,0,r") = h% la (r",0,r). (2.8) 

Thus, interchanging the summation variables a and a", and the variables r 
and r", we obtain 



n a e a n a „e a „ j d(r" - r) /fi a „(r, 0, r") 



fr — r") • r 



(r - r") 2 



= n a e a n all e a » f d(r" - r) /£ Q „(r, 0, r") — * 

a,a" [ r r ) 

= 2 ^ n a e a n a „e a n J d(r" - r) fo2, Q „(r, 0, r"), (2.9) 

a, a" 

where the last line is the half sum of the two first ones. 
Using (J21D in d2HD gives 

-2/3 _1 ^n Q / dr h aa >(r) = TTp b J2 n a"e a " / drr 2 Vo"( r ) 

a J a" J 

+e a <J2 n c* e a / dr h aa >(r) 

1 f 

+ - n a e a n a »e a/ : J dr d(r" - r) h^l, a „(r, 0, r"). (2.10) 

a, a" 

For the second term in the rhs of (12.101) . perfect screening [I] gives — e 2 /. For 
the last term in the rhs of (I2.10p . perfect screening gives 
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+ \ n a e a n a ne a n J dr d(r" - r) h®, a „(r, 0, r") 

a,a" 

= -- ^n a ne a n{e a > + e Q ») J dr h a > a "{r) 

a" 

= \ e a' ~ Y n a" e l" dr K> a "{r) . (2.H) 
La" 

Thus ( 12.101) becomes the general second-moment sum rule 

- (3Txp b ^n a e a f dr r 2 h a , a (r) = ]:^n a (A - (3e 2 a ) f dr h a > a (r) - \(3e 2 al . 

(2.12) 

By multiplying ( 12.12ft by n a /e a > and summing on a', one recovers the usual 
Stillinger-Lovett sum rule [3]. But (12.121) is a stronger sum rule. 



2.2 Guest charge in a one-component plasma 

We come to the case of a mixture of two species, with charge e% — e and 
density n 1 , charge e 2 = Ze and density n 2 , respectively. We choose a' = 2 
in (I2.12p . For dealing with one guest charge Ze only, we set n 2 = 0, rii = n, 
—pb = ne; n times the integral of h 2 i is —Z, by perfect screening. The sum 
rule ( 12.121) becomes 

(3irn 2 e 2 J dr r 2 h 2 i(r) = -2 

Since p(r\Ze,0) = neh 21 (r), (12TT5D is (F2TT]) . 

This result (12 .ip can also be retrieved by a different method in the next 
subsection. 



/ (3e 2 \ .pe 2 



(2.13) 



2.3 Another derivation 

(12.11) can be derived in another way if we assume that this second moment 
can be expanded in integer powers of Z. 

In the limit of small Z, the term linear in Z in (12. ip can be obtained by 
linear response theory. Indeed, if we introduce a guest charge Ze, located at 
the origin, into an OCP, the additional Hamiltonian is 

H' = Ze(j>(0), (2.14) 

where 0(0) is the microscopic electric potential created by the OCP at the 
origin. To first order in Z, the charge density at r is 

p(r\Ze, 0) = -/3(p(r)Ze0(O)) T = -Ze/3 J dr v(r')(p(r)p(r')} T , (2.15) 
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where (• • •) T denotes a truncated average. We define the Fourier transforms 

as 

/>) = jdr exp(ik-r)/(r). (2.16) 
Then, the Fourier transform of p(r\Ze,0) is 

p(k\Ze) = -pZe^S(k), (2.17) 

since the Fourier transform of v(r) is 2n/k 2 and the Fourier transform of the 
correlation of charge densities is S(k). For small k, S(k) has the expansion 



k 2 (1 - (3e 2 /A)k 
2^(3 4vr%/3 2 e 
Therefore, we get the zeroth moment 



Jdrp(r\Ze,0) = -Ze, (2.19) 

in agreement with equation (1.20) in [I], and the part linear in Z of the second 
moment (12. ip . 

It may be remarked that the k 4 term of (12.181) is related to the com- 
pressibility, which is exactly known only for the 2D OCP [5J. Therefore, no 
extension to 3D, with a closed result, seems possible. 

In the opposite case of large Z, the impurity expels the mobile parti- 
cles from a large region around it, leaving only the background. Essentially, 
p(r\Ze, 0) = —ne for r < R, where R is some large radius, and p(r\Ze, 0) = 
for r > R (there is a transition region [6] of width of the order n -1 / 2 , but 
in the limit of large Z, it gives a correction of lower order in Z). The ra- 
dius R is determined by the perfect screening condition (I2.19j) which gives 
R 2 = Z/{im). The second moment is 

j dr r 2 p(r\Ze, 0) = -nenR 4 /2 = -Z 2 ^—, Z -> oo, (2.20) 

which is the Z 2 term of (12. ip . and this is the highest-order power of Z in the 
second moment. 

The same argument extended to 3D gives 

/ dr i*PW*, 0) = -( 3 ^) 5/3 Z - °°" ( 2 - 21 ) 

Therefore, in 3D, the second moment is not a polynomial in Z, and no exact 
formula valid for any Z can be obtained by the present method. 
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2.4 Guest charge in a two- component plasma 



A sum rule for the TCP can also be obtained from (I2.12p . Now, we consider 
a mixture of three species, with charge e\ = e and density rt\ — n + , charge 
e 2 = — e and density n 2 — n_, charge e 3 = Ze and density n 3 , respectively. 
There is no background (p b = 0). We choose a' = 3 in (12.121) . Finally, for 
dealing with one guest charge Ze only, we set n 3 = 0, n + = n„ (neutrality); 
the system is stable against collapse if (3e 2 < 2 and (3Ze 2 < 2. We call 
n = n + + n_ the total density of the TCP. In f 1 2 . 1 2 j) appears the quantity 

n+h 31 (r) + n_h 32 (r) = n(r\Ze, 0) - n, (2.22) 

which is the excess density around the guest charge Ze. Then, f 1 2 . X 2 j) becomes 
a sum rule for the zeroth moment of this excess density: 

J dr [n(v\Ze, 0) - n] = Z 2 j^. (2.23) 

This result is a generalization of the compressibility sum rule [7] 

/On 
dr[n(r|±e,0)-n] = ^-l (2.24) 

with the use of the exact equation of state j3p = n{l — Pe 2 /4), where p is the 
pressure. 



2.5 Mixture without a background 

In the case p b = 0, another derivation of f)2.12p is possible starting from the 
known equation of state [8] 

The M-component plasma may be described in the grand-canonical ensemble, 
with M chemical potentials p a (actually [9], the system turns out to be neutral 
in the thermodynamic limit, and M — 1 chemical potentials would suffice for 
determining the state of the system; but here it is more convenient to use M 
chemical potentials). The pressure p is given by (3p = lim(l/K) InS, where 
V is the volume (here area) of the system, S is the grand partition function, 
and lim is the thermodynamic limit. Taking the partial derivative of (I2.25P 
with respect to (3p a i gives 




(2.26) 



which is f l2~T2D with p b = 0. 
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3 Guest charge and potential fluctuations 



Putting a guest particle of charge Ze at the origin r = 0, the original Hamil- 
tonian Hq of the infinite Coulomb system modifies to H = Hq + Ze<p(0), 
where 0(0) is the microscopic electric potential created at the origin by the 
Coulomb system. The charge density around the guest charge, at point r, is 
thus expressible as 

(p(r)exp[-/^e0(O)]) 
p[T\Ze,0) = — - r — ~ n1 , , (3.1) 



\ CX P 



-f3Ze<p(0) ) 



where (• • •) denotes the thermal average over the homogeneous system with 
the Hamiltonian Hq. 

Let f/fze denotes the excess (i.e., over ideal) chemical potential of the guest 
charge, i.e. the reversible work which has to be done to bring the guest particle 
of charge Ze from infinity into the bulk interior of the considered Coulomb 
plasma. By the coupling parameter technique [7] , this chemical potential can 
be represented in terms of the charge density (13.11) as follows 

/4 x e = eJ^dZ'J dr v(r)p(r\Z'e,0). (3.2) 

With regard to the representation (13.11) . pf e can be expressed as 

-K=f'\(^». (3.3) 
J o (exp(x0); 

Here, since the thermal averages are point-independent, we use the notation 
= 0(0). ' ' 

Let us recall some basic information about the cumulant expansion. Let 
be a random variable with the probability distribution P(4>). The cumulant 
expansion is defined by 

/ 00 x l - \ 

(exp(x0))=exp(£-(0') c J, (3.4) 

where x is any complex number and (4> l ) c are the cumulants. They are com- 
binations of the standard moments (0 Z ). Differentiating the equality (13. 4p 
with respect to x gives 

J OO l ( OO I \ 1 oo I 

Sg|^-xp(g|<^)^E|^), (3.5) 

Equating the coefficients of the same power of x in both sides of (13. 5p gives 
the recursion formula 
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The first cumulants read 



(4>)c = (0), 

(0 2 > c = <0 2 >-<0> 2 , (3.7) 
(0 3 ) c = #)-3(0 2 )(^) + 2(M 

etc. In the theory of fluids, the cumulants of type (13. 7p are referred to as 
truncations, and therefore we shall use the notation {4> l ) c = (4> l ) T - 
Since it holds 

<^EWl> = iL,n(e X pW)), (3.8) 
(exp(x0)) dx 

the excess chemical potential (13. 3p is expressible as 

-/3/i| x e = ln(exp(-/5Ze0)). (3.9) 

Based on the recapitulation in the above paragraph, /i^ e is expressible either 
in the form of a cumulant expansion 

- = E ^P^r, (3-10) 

or in the form of the standard moment expansion 

exp = (exp(-/3Ze<j>)) = 1 + £ ^-^)- (3.11) 

It stands to reason that the expansions (13.101) and ( 13. lip are valid provided 
all moments exist. We conclude that the knowledge of the excess chemical 
potential of the guest particle with an arbitrary charge provides the exact 
information about all moments of the electrostatic potential at a point of the 
infinite homogeneous Coulomb system. 

Going to the infinite system via the thermodynamic limit of a finite sys- 
tem with a disc geometry [TO], the fluctuations of the potential at any point 
become infinite due to the presence of dipoles near the boundary. Here, the 
potential moments are defined directly for an infinite space, without the pres- 
ence of a boundary. This corresponds to going to the infinite system via the 
thermodynamic limit of a finite system, e.g., with periodic boundary condi- 
tions, formulated on the surface of a sphere and so on. We thus expect that 
the average potential at a point is equal to zero and all its moments are finite. 

Since in 2D the potential (11.21) is dimensionless, <p has the dimension of 
the elementary charge e. It is therefore useful to introduce the dimensionless 
microscopic quantity ip = <p/e with the probability distribution P{ip). Setting 
in (13.111) [3Ze 2 = ik, one gets 

exp (-/3/4 x J | ^ = (exp(-i^)) = / d^ e-^P(^) = P(k), (3.12) 



11 



where P(k) is the Fourier component of the -^-distribution. The original 
probability distribution P(ip) can be obtained by the Fourier inversion of this 
relation 

/OO Ah 
oo -e^e W {-^ e )\, Ze2=lk . (3.13) 

All that has been said in this section is valid also for v being the pure 
Coulomb potential plus any type of short-distance regularization. 



4 High-temperature limit 

The high-temperature (weak-coupling) limit of Coulomb systems is described 
rigorously by the Debye-Huckel theory [HI [12]. In 2D, the two-body Ursell 
functions U of charged species are given by [T3] 

U aa >(r,r') = n^,(r,r') - n a n a , = -e a n a e a ,n a > [3K (k\t - r'|), (4.1) 

where K is a modified Bessel function [H] and k = (2ttP J2 a e a n aY^ 2 is the 
inverse Debye length. 

The potential-potential correlation function can be calculated directly 
from the definition 

(0(O)0(r)) T = J drx v(r - n) J dr 2 t;(r 2 )(p(r 1 )p(r 2 )) T 

= J drx v(r - n) J dr 2 u(n - r 2 )(p(0)/3(r 2 )) T . (4.2) 
Using for the Coulomb potential the expansion in polar coordinates 



oo 1 / r \i 



v{T l -v 2 ) = -\n\r l -r 2 \ = -\nr > +Y,l\ — \ cosZ(0i-0 2 ) (4-3) 



i=i Z V r >, 



with r < = min{r 1; r 2 } and r > = max{r!,r 2 }, and taking into account the 
screening sum rule [3] 

J dr 2 (p(0)p(r 2 )) T = 0, (4.4) 

the second integral on the rhs of f!4.2p can be expressed as 

J dr 2 u(n - r 2 ) (p(0)p(r 2 )) T = - jH dr 2 2vrr 2 In (p(0)p(r 2 )) T (4.5) 

Considering the charge correlation function 

(p(0)p(r 2 )) T = ^ e a e a > \u^l,{r 2 ) + n a S aa 'S(r 2 ) 



a, a 



^o(«r 2 ) + — -<f(r 2 ) (4.6) 



(2tt) 2 /? uv zy 2tt/3 
12 



in equation (14.51) implies, after an integration by parts, 

2 

J dr 2 viv, - r 2 )(p(0)p(r 2 )) T = JL-K q {kt x ). (4.7) 

Inserting this relation into (14. 2[) and applying once more the expansion H4.3[) 
results into 

/3(0(O)0(r)) T = -lnr-K {Kr). (4.8) 

This procedure will be repeated, without going into details, also in the cases 
treated in the next sections. 

The result (14.81) has the correct large-distance asymptotic [15] 

/3(0(O)0(r)) T r:roo -lnr. (4.9) 

In the zero-distance limit r — > 0, using the expansion 

K (x) = -C-\n(x/2) + 0(x 2 \nx) (4.10) 

with C being the Euler number, the one-point second-moment fluctuation 
formula for the potential reads 

P(4>y = C + \n(K/2). (4.11) 

One can obtain the last result in an alternative way by considering the 
charge density induced around the guest charge [1] 

p(r\Ze,0) = -Ze^K ( K r). (4.12) 
Then, according to (13. 21) . 

-Pfif e = -(3e 2 dZ' Z'— dr 2nr In rK (Kr) 
Jo 2n Jo 

= ^^[C + Hk/2)}. (4.13) 
With regard to the cumulant expansion (13.101) . we recover the previous result 

(BHD- 

From ( 13.101) and ( 14.131) . all the higher-order truncated moments ((p l ) T with 
/ > 3 vanish in the Debye-Huckel limit; this indicates a Gaussian distribution 
for the one-point potential in this limit. We shall return to this problem and 
present all truncated potential moments, for the TCP, in a high-temperature 
limit going beyond the the Debye-Huckel limit, in Sect. 6. 
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5 One-component plasma at (3e 2 = 2 

The 2D OCP is exactly solvable in terms of free-fermions when the dimen- 
sionless coupling constant (3e 2 has the special value 2 p^6l fTT] . In the thermo- 
dynamic limit, the two-body Ursell function of mobile particles at distance r 
is 

U(r) = —n 2 exp 7mr 2 ) , (5-1) 

where n is the particle density. All many-body Ursell functions are known at 
the free-fermion point, too. 

The potential-potential correlation function can be calculated in close 
analogy with the previous steps outlined between Eqs. (14.21) - (14. 8p . Sub- 
stituting the charge correlation function 

(p(0)p(r 2 )) T = -eV exp(-7mr 2 ) + n5(r 2 ) (5.2) 

into the relation H4.5[) and using an integration by parts, one gets 

/* €. 2 Tl 

J dr 2 v{r x -r 2 )(p(0)p(r 2 )) T = _r(0, vmr 2 ), (5.3) 
where _ 

poo 

„x— 1„— s 



T(x,t) = j ds s^e- 8 (5.4) 
is the incomplete Gamma function. From (14.21) . one thus obtains 



/3(0(O)0(r)) T = - In r + - [e^ 2 - (1 + vmr 2 )r(0, vmr 2 )J . (5.5) 

This result has the correct large-distance asymptotic (14. 9j) . In the zero- 
distance limit, it yields 

P(iY = i[l + C + ln(7m)]. (5.6) 

Note that the large-distance behavior H4.9I) is universal, while the zero-distance 
limit (14.111) or (15.61) depends on the coupling constant /3e 2 . 

All potential moments are available for the present system due to the 
knowledge of the induced charge density around the guest charge [TH [1] : 

p(r\Ze,0) = -en \ Z > 0. (5.7) 

By using the relation (13.21) . one obtains after some algebra [IB] 

- m* = %• t 1 + ln (™)] - I" dZ ' Z '^ 1 + (5-8) 

2 Jo 
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where ip is the psi-function defined by 



Mx) = ^-hxVix). (5.9) 
ax 

Its Taylor expansion around x = 1 reads [H] 

oo 

+ x) = -c + ^(-i)'c(O^ 1 , (s.io) 

Z=2 

where 

OO 1 

C(0 = Eu (s- 11 ) 

fe=l K 

is the Riemann zeta function. Considering the expansion f)5.10p in (15.81) gives 

- me = T [1 + ^ + !*(*»)] + £ { -^—C(l - !)• (5.12) 
The comparison of this expansion with the cumulant expansion ( 13.101) implies 

(0 2 ) T = j[l + C + Hnn)}, (5.13) 

(0') T = ^(/-1)!C(/-1), ^>3. (5.14) 

Note that the second-moment formula (I5.13P is identical to the previous one 
(15.61) derived by the direct calculation from the definition. 



6 Two-component plasma 

6.1 Collapse point f3e 2 = 2 

The 2D TCP of ±e charges is mappable for the special value of the coupling 
constant f3e 2 = 2 onto the Thirring model at the free-fermion point [19], [20] . 
Although this coupling corresponds to the collapse threshold for the pointlike 
particles, and therefore for a fixed fugacity z the particle density is infinite, 
the Ursell functions are well defined. Their two-body forms read 

U±,±(r) = ~ (^j K*(mr), U±, T (r) = (^-) tf*(mr), (6.1) 

where m = 2ttz. All many-body Ursell functions are also known. 
Substituting the charge correlation function 

(p(0)p(r 2 )) T = -2e 2 — \Kl{mr 2 ) + Kl{mr 2 )\ (6.2) 
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into the relation (14. 5p and integrating by parts leads to 

dr 2 u(n - r 2 )(p(0)p(r 2 )) T = e 2 — K 2 (mn). (6.3) 
From (14.21) . one finds that 



[mr 



\2 



/3(0(O)0(r)) T = -lnr + 2fsT 2 ( 

(6.4) 

This result has the correct large-distance asymptotic (j4.9j) . In the zero- 
distance limit, it gives 

/3(0 2 ) T = l + C + ln(7rz). (6.5) 
6.2 Stability region < (3e 2 < 2 

The system of pointlike ±e charged particles is stable against the collapse of 
positive-negative pairs of charges provided that the corresponding Boltzmann 
weight exp[/3e 2 v(r)] = r~^ e can be integrated at short 2D distances, i.e. 
when f3e 2 < 2. The equilibrium statistical mechanics of the neutral TCP is 
usually studied in the grand canonical ensemble, characterized by the particle 
fugacities z + = z_ = z. The full thermodynamics of this system is known 



In the stability range of (3e 2 < 2, the grand partition function H(z) of the 
2D TCP can be turned via the Hubbard- Stratonovich transformation (see, 
e.g., Ref. [23]) into 

\ ST>(pex$[-S{z)\ 



" v ' JVipexp[-S{0)Y 

where 



S(z) = J dr 



-(Vy?) 2 - 2zcos(b<f) 



(6.6) 



(6.7) 



167T 

is the Euclidean action of the (1 + l)-dimensional sine-Gordon model. Here, 
(p(r) is a real scalar field and / T>(p denotes the functional integration over this 
field. The sine-Gordon coupling constant b depends on the Coulomb coupling 
constant via 

»=)/?■ < 6 ' 8 > 

The fugacity z is renormalized by the diverging self-energy term exp[j3v (0)/2] 
which disappears from statistical relations under the conformal short- distance 
normalization of the exponential fields [2TT [22] 

^ e iV(r) e -i^(r')^ G _ | r _ r /|-46 2 ag | r _ r /| ^ Qj 
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where (• • -) s g denotes the average with the sine-Gordon action (16.71) . The 
species densities are expressible in the sine-Gordon format as follows 



n± 



sG- 



71- 



(6.10) 

n/2 is ensured by the obvious 



The charge neutrality of the system n. 
symmetry relation (e lbip ) S G = (e~ l6</3 ) S G- 

The excess chemical potential of the particle species forming the plasma 
is given by 



exp( 



"/Vie) = ^ 



(e ±ibv ) 



sG- 



(6.11) 



It was shown in Ref. [24] that the excess chemical potential of a guest charge 
Ze immersed in the plasma is expressible in the sine-Gordon format as follows 



exp(-/?/4 x e ) = (e^) sG 



(6.12) 



When Z = ±1, one recovers the previous result (16.111) valid for the plasma 
constituents. Due to the symmetry relation (e m{p ) S G = (e _ia ^%G valid for any 
real-valued a, it holds that [f^e = ^-Ze- 

The (1+1) -dimensional sine-Gordon model is an integrable field theory 
|25j . Due to a recent progress in the method of the Thermodynamic Bethe 
ansatz, a general formula for the expectation value of the exponential field 
(e ia< ^) was derived by Lukyanov and Zamolodchikov [26]. In the notation of 
equation (I6.12p . a = Zb, their formula reads 



JZbip 



)s 



G 



with 



h{Z) 



00 dt 

T 



nzT{l - b 2 

T\b 2 



smh 2 (2Zb 2 t) 



2 ,- l (Zbf/(l-b 2 ) 



exp [h{Z)\ 



(6.13) 



2 sinh(6 2 t) sinh(t) cosh[(l - b 2 )t] 



- 2Z 2 b 2 e 



(6.14) 



The interaction Boltzmann factor of the guest charge Ze with an opposite 
plasma counterion at distance r, r~@ e is integrable at small 2D distances 
r if /3\Z\e 2 < 2, i.e. \Z\ < 1/ (2b 2 ); this is indeed the condition for the integral 
(I6.14p to be finite, so that the couple of Eqs. (I6.13P and (16.141) passes the 
collapse test. Finally, using eqs. (16.131) and (16.141) in (I6.12p . one arrives at 



Pi* 



ex 
Ze 



b 2 



b 2 



In 



vr^rfi 



T(b 2 



h{Z). 



(6.15) 



We have to keep in mind that b 2 = j3e 2 /4. 
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Comparing the cumulant expansion (13.101) with the result (I6.15P , in which 
the integral h(Z) (16.141) is expanded in powers of Z, one gets the explicit 
forms of the potential moments: 



-2 



V = - In 

86 2 (1 - b 2 ) 

e 2 f°° dt 
+ — 



, nzT(l - b 2 ) 



Y(b 2 



4 Jo t 



t 2 I 



sinh(6 2 t) sinh(t) cosh[(l - b 2 )t] b 2 



(6.16) 



{<p2l)T ~ 4 Jo dt sinh(6 2 t)sinh(t)cosh[(l-6 2 )t]' 1 2 ' 3 ' ' ' ' ' (6 ' 17) 

The odd potential moments vanish for the symmetric TCP. 

In the high-temperature limit (3e 2 — > (b 2 — > 0), (16.151) taken with z ~ n/2 
reduces to the previous one (14.131) ; one retrieves the second moment (14. lip 
and that all higher moments vanish, as it should be. From (I6.17p . in the limit 
b 2 — > 0, one finds 

P(4> 2l f = e 2 ^8^^(2l - 2)!C(2/ - 1), 1 = 2,3,.... (6.18) 

These expressions go beyond the Debye-Htickel limit of (16.151) . 

At the collapse point (3e 2 = 2 (b 2 = 1/2), the second-moment formula 
(I6.16P reproduces the previous result (16. 5p and the higher-order moments 
(I6.17P take forms 

(^) T = e2/ |( 2/ " 1)'C(2Z - 1), I = 2, 3, . . . . (6.19) 

All potential moments are finite also in the collapse region, up to the 
Kosterlitz-Thouless critical point j3e 2 = 4 {b 2 = 1). We conjecture that, in 
the case of the hard-core regularization of the Coulomb potential, the obtained 
result correspond to the limit of a vanishing hard core. 

We end up this section by a comment about the possibility of a relationship 
between the electrostatic potential <fi and the sine-Gordon field variable ip. 
This relationship was suggested in many articles, see, e.g., Ref. [27]. The 
comparison of Eqs. (13.111) and (I6.12p implies 

(vAg = (-l) l m l (t 21 )- (6.20) 

This means that, in view of one-point fluctuations, the fields and <p differ 
from one another only by an irrelevant scaling factor. On the other hand, 
the large- distance asymptotic of the potential-potential correlations (14.91) is 
fundamentally different from the one of {ip(0)ip(r)) T The latter two-point 
correlation function has, like in every massive field theory, a short-range ex- 
ponential decay as r — * oo. We conclude that the electrostatic-potential 
interpretation of the sine-Gordon field is not correct. 
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7 Conclusion 



The general study of a mixture of M species of mobile particles, which may be 
embedded in a uniform background, is simpler in two dimensions; the BGY 
hierarchy suffices for deriving the general sum rule (12.121) relating the second 
moments and the zeroth moments of the two-body correlations. Further work 
should be possible about this mixture. 
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